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ABSTRACT 



Rational solutions of the Painleve IV equation are constructed in the setting of 
pseudo-differential Lax formalism describing AKNS hierarchy subject to the addi- 
tional non-isospectral Virasoro symmetry constraint. Convenient Wronskian repre- 
sentations for rational solutions are obtained by successive actions of the Darboux- 
Backlund transformations. 

1 Introduction 

Two main observations are put together in this paper with the purpose of developing 
a systematic method rooted in Darboux-Backlund techniques within the pseudo- 
differential Lax formalism to generate rational solutions of the Painleve IV equation. 
First, we observe that the AKNS hierarchy subject to the additional non-isospectral 
Virasoro symmetry constraint (also known as the string equation) reduces to the 
Painleve IV equation. Secondly, we make use of the fact that the Darboux-Backlund 
(DB) transformations commute with the additional-symmetry Virasoro flows [2] and 
thus these transformations can effectively be used to construct all known rational 
solutions O O [TOl [12] of the Painleve IV equation (14.31) out of a few basic ones. 



2 AKNS hierarchy with string condition 

To introduce the AKNS Lax formalism with additional Virasoro symmetry flows we 
follow |2j and define the pseudo-differential Lax hierarchy as described by a dressing 
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formula 



The associated t„-flows of this hierarchy: 



dr 



Ll(r) 



reproduce for n = 2 the AKNS equations 



d d 2 

TT-q + qxx-'^q r = 0, — r - r^x + '^qr = . 



[2.1] 



(2.2) 



(2.3) 



L2 

Furthermore, it is useful to introduce as in [3J a squared eigenfunction p such that 

Px = -2rq 
and, which in particular, satisfies |3] : 

d 



dt 



-p = 2 {-qrx + rqx 



(2.4) 



The above formalism is augmented by the Orlov-Schulman operator [T 



M=W (^Itidl-^^ W-^ 



conjugated to the Lax operator according to the commutation relation [L , M] = 1. 
The additional Virasoro symmetry flows are then defined as 



dk,nL — — 



n—l T k 



(2.5) 



and they commute with the usual isospectral flows as defined by dL / dtn = [L^ , L] 
or (Q. 

For n = 1 formula 02.51) yields 



{dk,iL)_= [{ML'')^,l\ +{L') 
which can be rewritten as 



(2.6) 



fc-i 



{dk,iL)_ = {ML'')^ {r)d-^q-rd-^ {ML'^)*^ {q) + J2L'''''\r)d-^ {L*y (q) . (2.7) 

For k = 0, 1, 2 these flows form the s/(2) subalgebra of the Virasoro algebra and 
preserve the form of the Lax operator of the AKNS hierarchy. Subsequently, they 
can equivalently be described by their action on the eigenfunction r and the adjoint 
eigenfunction q : 

{drL)_ = (drr) d-^q + rd-^ (9,g) (2.8) 
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with dr = dk^i {k = 0, 1, 2). Explicitly, one finds [2] : 

Boat = (M)^ (r), do,iq = - (M); (g), (2.9) 

ai,ir = (ML)^ (r) + (1 - z/)r, Bi^iq = - iML)l (q) + uq, (2.10) 
d2,ir = {ML')^ (r) + L(r), Ss.ig = - (ML^); (g) + L*(g). (2.11) 

Note arbitrariness in equation (12.10p expressed by a parameter i/. 
In what follows we impose the symmetry conditions: 

d,,,r = 0, d,,,q = 0. (2.12) 

Since 

ML = W (^Itidl^j W~\ (ML)* = W'^* {-d^x + dl2t2 + ...)W* 
= -1 + W^'* {-xd, + 2t2dl + ..)W* 

it follows that 

(ML)^ = xd^ + 2t2 {L)l + . . ., - {ML)l = l + xd^- 21^ {L*)l + ... 
Thus in view of equation (12.101) the condition (12.121) amounts to 

d d 
-xq^ - 2t2Tr-q = q + J^q, xr^ + 2t2Tr-r = -r + ur 

dt2 

The model we consider is obtained by applying the string equation (12.131) to eliminate 
t2-dependence from r, q, p and setting the higher flows to zero. Setting t2 = —1/4 and 
eliminating t2-6^ow dependence from eq. (12.131) by inserting values from (I2.3l) - (l2.4p 
yields : 

-xqx + ^ {-qxx + 2gV) = q + uq, xr^ - ^ (r^.^ - 2gr^) = -r + vr 
p + xpn, = p- 2xrq = q^r - qr^, 

Multiplying the first equation by and second by g^. and adding them one obtains: 

{qxrx - q^r^)^ = '^{qxT - qr^) - 2v{qr)^ 

or after integration : 

qxfx = q r + 2xp — 2iJrq — [p — ^ ) , (2.15) 

where we have set the integration constant to /i^ — z/^, with p being an extra parame- 
ter. In what follows p will, in addition to u, parametrize solutions of the Painleve IV 
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equation. The origin of /i is very transparent in an algebraic approach to integrable 
models subject to the scaling condition [1]. 

Dividing the first of eqs. (12.141) by q and second by r and summing them yields: 

2. = xf-^ + ^V-f^-2^V2gr = 



q r J 2 \ rq rq 

x[2x-P-)-'-(^^-2'-^)+2qr, 
rq J 2 \ rq rq J 

where in the last line we used the third of eqs. (I2.14p . Inserting q^Vx from eq. (I2.15P 
produces after multiplication by rq : 

2x^rq + xp — -^{rq)xx — 4:h'rq + 3(rg)^ = p^ — v"^ . 

Inserting in the above rq = —px/2 results in an equation entirely expressed in terms 
of only one variable p : 

1 3 
- x'^px + xp + -pxxx + 2iypx + -pi = p^ -iJ^ . (2.16) 

One recognizes in the above equation the special case of the Chazy I equation [7]. 
After multiplication of eq. (I2.16P by pxx we can rewrite the resulting equation as a 
total derivative, which after an integration becomes 

pL = 4 {xpx - pf - 2pl - 8z/p2 + 8(/i2 - u^)px - 8C , (2.17) 

with C being an unknown integration constant. Another well-known form of this 
equation is [9] : 

3 

pIx = 4 {xpx - p)^ - 2 JJ {px + Vi) . 

i=l 

By comparison with eq. (I2.17P one obtains 

vi + V2 + V3 = Au, V1V2 + Viv^i + V2V3 = -A{p'^ - ly'^) , V1V2V3 = AC . 

If we set the integration constant C to zero then eq. (I2.17P can be simplified as 
follows : 

pL = 4 (xpx - pf - 2p, [p, - 2(p - u)] [p, + 2(p + u)] (2.18) 
or in a slightly different and more convenient form : 

(2 (xp, - p) + Pxx) (2 (xp, - p) - Pxx) = 2px [px - 2(p - z/)] [p, + 2(p + z/)] . (2.19) 

First, one notes that equation (12.191) is manifestly invariant under p — > — p and thus 
a solution of equation (12.191) with one value of p solves this equation for — p as well. 
Next, one notes that the left hand side of equation (I2.19P remains invariant under 
substitution p = p+kx with k being a constant. For k = 2{p — v) and k = —2{p + i>) 
the right hand side can be given the form 2p,[p, — 2(p — ^')][p, + 2(p + u)] with 
u = 3p/2 — Z//2, p = ±(p + z/)/2 and u = — 3p/2 — z//2, p = ±(p — z/)/2, respectively. 
Hence for these values of k the transformation p ^ p takes the "old" solution of 
equation (12.190 to the "new" solution of equation (I2.19P with the new parameters 
p, u. 
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3 Basic Polynomial Solutions 



One finds by inspection that there exists a class of polynomial solutions of equation 
(I2.19P for p being of the form p = bx^ + ex. With both coefficients, b and c, being 
non-zero the polynomial solution is : 

8 4 1 
p= — x^±-x, 11^ = -, z/ = =1=1 . (3.1) 
p 27 3 ' 9' ^ V ; 

For c = 0, one finds that the unique non-zero solution requires b = 8/27 with 
zy = 0, /i2 = 4/9 : 

P=^^', P'=g, ^^ = 0. (3.2) 

Finally, one notes that setting 6 = causes both xpx — p and pxx to vanish together 
with the left hand side of eq. (I2.19p . Thus p = cx is a solution to eq. (12.191) for 
three values of c, namely c = 0,2{p — u), — 2(/i + u) for which the right hand side of 
eq. (12.191) vanishes as well. 

In a text below, the polynomial solutions of equation (I2.19P will serve as seeds 
of chains of the Darboux-Backlund transformations. 

For p = cx the product rq is a constant, — c/2. In addition, qxr—qr^ = p—2xrq = 
2cx and thus rqx = cx, r^.g = —cx. Multiplying eq. (12.150 by r and then by q one 
gets from rqx = cx, r^g = —cx that 

cxTx = (2cx^ + c^/4 + cu — (/i^ — z/^)) r, —cxqx = (2cx^ + c^/4 + cu — (/i^ — z/^)) q 

For c = 2{p — u) , —2{p + u) the above equations simphfy to 

xrx = 2x^r, —xqx = 2x^q , — > r = kic^^ , q = k^ie~^^ 

and as long c = —2kik^i ^ the above solutions do not lead to non-zero solutions 
of the Painleve IV equation. 

We now consider the case of 6 = c = 0. Then rq = with either r = or g = 
and it follows from (I2.15P that p^ = z/^. 

Another useful expression for eqs. (I2.14p is : 

(^e'^'V^j + 2(z/ - l)e~''V = 2(rg)e"'''V, (e^'^x^ + 2 (z/ + 1) e^^'g = 2{rq)e^^ q 

(3.3) 

For g = and v — 1 = n being a positive (or zero) integer the first equation 
of eqs. (13. 3p is Hermite's equation, alternatively known as r^^ — 2xrx + 2nr = 
with r{x) = Hy_i{x) being a solution. Thus one obtains, r = Hq = 1 loi v = 1, 
r = Hi = 2x for u = 2 (the solution we have already studied) and r = H2 = 4x^ — 2 
for z/ = 3 etc. 

For — z/ > being a positive integer a substitution 

r{x) = e^^ g{x) 

leads to the following equation 

gxx + 2xgx + 2vg = gxx + 2xgx - 2mg = 0, m = -u = 0,1,2, .. . 
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which has a polynomial solution 



v2 



gm{x) = H^x) = e-^ = (-i)"ff^(ix) (3.4) 

with gm=o = 1, gm=i = 2x, gm=2 = 2(1 + 2x^) and gm=3 = (3 + 2x^)4x etc. To 
summarize : 

.(.jJ^i. 'f"*"'^-;'^- (3.5) 

I e H_i,[x), II z/ < U IS mteger 

for g = 0. Substituting 

g = e-^V(a;), 

into (13.31) with r = one gets the following equation for / : 

f.. - 2x/, + 2uf = 0, 

which is a Hermite's equation for Hn{x) for u = n > 0. Thus f{x) = Hy{x) = Hn{x) 
with n = 0, 1,2, . . . and q{x) = exp{—x^)H^{x) as long as z/ is a positive integer. 
For z/ + 1 = — m, where m = 0, 1, 2, . . . the equation for q, 

Qxx + 2xqx - 2mq = 

is identical to equation for g written above and has the same solutions. To summa- 
rize: 

J e-'^'if^, if 1/ > is integer 
g(x) = < ^ (3.6) 
\H_^_i{x), if z/ < is integer 

for r = 0. 



4 Painleve IV equation and its Hamiltonian for- 
malism 



In terms of variables: 



y = -—-2x, w = — -2x, z = -2rq + 2(a + u) (4.1) 
q r 



we can express equations (12.141) as 



Ux = z + y + 2xy — 2fi, Wx = —z — w — 2xw + 2fi 



Z^ X „ -2 



Zx = 2rq{y — w) = h wz — 2{jj, + u)w + 2/i— 



h wz — z,\^ -t- z/ru; -t- ZyU— /a r)\ 

2w w \^-^) 

z^ z 

= TT ^ + + - 2/^- 

2l/ y 

Applying one more derivative on the first equation and eliminating z yields the 
Painleve IV equation: 

13 

= i^vl + 77?/' + W + 2(^2 + zy + l)y - 2^ (4.3) 
2?/ 2 y 

6 



Repeating the same procedure for the equation with w yields almost identical equa- 
tion (note however a shift in u as compared to eq. (14.31) ): 

13 

Wxx = — wl + -w'' + Axw^ + 2{x^ + v-l)w- 2— (4.4) 
2w 2 w 

As found by Okamoto [13j, equations (14.21) for a pair of variables {y,z) possess a 
Hamiltonian representation. Define namely 

H = -2P^Q + (Q2 + 2sQ - 2/i) P + ^ (/i + z/)Q . (4.5) 

Then the Hamilton equations are: 

OH 

= T77T = -4gP + + 2xQ - 2/i (4.6) 

riM 1 
" ~ 9g " " - 2xP - - (/i + z/) . (4.7) 

They agree with eqs. (14. 2 p for Q = y and z = —4QP or P = —z/Ay. Due to the 
Hamilton equations it holds that 

A.H = = 2QP = -^z = rq-{u + i^) = -^p,. - (i/ + /x) . (4.8) 
Thus, up to a constant, 

H = -^p-x{iy + fi) . (4.9) 

Plugging the product QP = — {p^ + 2{p + u)) / 4 into expression (14. 5 p for H and 
using eq. (14. 9 p one gets 

P {-rq - p + z/) + Q^rq + xrq = -^p 



or 



Define: 



Then 
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-P (p:,. - 2(p - z/)) - -Qp^ = - {xp^ - p) . 

P = P^ ^7 V Q = ~^Qo-. (4.10) 

p^ - 2(/i - z/) 



^oQo = -^Px {px - 2(p - v)) {px + 2(p + z/)) = ^ (4 (xp^. - p)^ - pI^) 

= ^ (2 (^^P^ - P) + Pxx) (2 (xp^. - p) - Pxx) 
d4 

where use was made of eq. (12.181) . Solutions to the above equation are easily found 
to be ^ ^ 

Po = - (2 {xpx - p) ± Pxx) , Qq = -{2 {xpx - p) =F Pxx) , 



which gives two answers for expression for the Painleve IV solution y = Q in terms 
of the solution to the p-equation fl2.19p 

y± = Q = ^{2 {xp, -p)± p,.,.) (4.11) 

The above ambiguity in signs can be explained by comparing with the third equation 
in (14.21) . Consider namely the difference ?/+ — ?/_ of two solutions from eq. (14.111) : 

Pxx ^x 

y+-y- = = 7^ = y - w 

Px 2rq 

where we used the third of equations in (14. 2p . Thus ?/+ = y, y- = w. It is natural 
to connect these associations by a Backlund transformation by defining a "barred" 
system such that 

w = = y+ = y = {2 {xp^ - p) + pxx) 
2px 

and thus the new p variable needs to satisfy 

P ^ Pxx P Pxx 

px 2p^ p^ 2px ' 

From the last of equations (12.141) one can obtain expressions for the right and the 
left hand side of the above relation in terms of q and r. This results in a condition : 

K 

(Inr)a, = -(Ing)^, -> g = — , K = const , (4.12) 

r 

which also easily follows from the equality w = y. The above relation resembles an 
effect of the Darboux-Backlund transformation to be described in the next section. 
Note that the above transformation maps w to y and therefore comparing with 
values of z/ in eqs. (14.31) and (14. 4p we conclude that it lowers u hj 2. 

For completeness, let us mention that as follows from eqs. (14. 7p ioi Q = y and 
P = —z/Ay, the function Y = —2P = z/2y also satisfies the Painleve IV equation 
(14.30 with the parameters u = 3p/2 — i//2 — 2, p = ±(p + z/)/2. By comparing with 
definitions (I4.10p and definitions of Pq, Qq we find that the transformation y —>■ Y 
is a superposition of the transformation p ^ p discussed below equation (12.190 with 
the Darboux-Backlund transformation y —>■ w lowering u by 2. In what follows we 
will only discuss the latter transformation together with its square-root. 



5 Darboux-Backlund transformations 

The Darboux-Backlund transformation is here introduced as a similarity transfor- 
mation with an operator T = rd^r^^ acting on the pseudo-differential AKNS Lax 
operator through : 

L = d^-rd-^q L = T LT-^ = d.^-fd-^q 
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A simple calculation yields 

f = r(lnr)^^ - r^g, q = -- . (5.1) 

r 

in agreement with relation (14.121) for K = —1. Taking a product of f, q one obtains 

= -2fq = px + 2{\iar)^^ . 
Thus the DB transformation (15.11) yields 

p = p + 2(ln r)x = p + 2w + 4:X = p + 2y^ + Ax 
= p + 2^ (2 (xp, - p) - p,.J + 4x = p + 2-^ + ^ + 2x (5.2) 

2px Px Px 

or including a reference to the (p, v) parameters 

(M,-2) ^ (M,.) + + ^ + 2x . (5.3) 

Px Px 

It is equally easy to formulate the adjoint Darboux-Backlund transformation which 
increases v by +2. The adjoint Darboux-Backlund transformation involves acting 
with 5* = qdxQ^^ on the pseudo-differential Lax operator through : 

L = dx-rd-'q L = S*-'LS* = {q-'d~\)L{q-'dxq)=dx-7d;'q 

with 

q = -q{\ia.q)xx + q^r, r = 

q 



q = -q(\ia.q)xx + q^r, r = - . (5.4) 



Taking a product of r, g one obtains 

Px = -2rq = p^ + 2{lnq)^^ 
Thus the DB transformation (15.41) yields 

p = p + 2(ln g)^ = p - 2?/ - 4x = p - 2t/+ - 4x 
= p - 2^ (2 {xpx - p) + Pxx) -Ax = p-2-^ + ^-2x (5.5) 

2px Px Px 



or 

i,,u+2) ^ M _ + Pff_^ _ . (5.6) 

Px Px 

Transformations and " when acting on variables J, J such that: 
J = -rq = p^/2, J = {\nq)x = -y -2x 
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take the following form [T] : 

G{J) = J+{\n{j+J^))^ G{J) = J + J, (5.7) 

G-\J) ^J- (In J)^ G-\J) ^ J + (In J)^^ - J,, (5.8) 

where we found it convenient to rewrite actions of ~ and " as transformations G 
and G-\ It follows from dEZD-O that: 

G{y) = y-{\n {y, + y^ + 2xy + 2v + 4))^ , G-\y) = y+{\n {y, - y^ - 2xy - 2u))^ 

(5.9) 

The DB transformations G^^ are equal to Murata's transformations first obtained 
in [H]. This is established after making use of the Painleve equation (14. 3 p in (15. 9p 
and identifying parameters 6, a from [llj with /i, —u — 1. 

There also exists a set of variables j, J related to J, J via a Miura transformation 
m such that 

J = ~3-j + - ; J = JJ (5.10) 
J 

In terms of variables j,j one can define a discrete symmetry transformation 

9(3) =3-— 9 {3)= 3 (5.11) 

3 

g-\j)^j+h g-\j)^j (5.12) 

3 

such that [T] 

^?' = G, (5.13) 

when applied on J, J defined by equation (15.101) . 

From the above relations one obtains simple transformation rules: 

g{y) = y-{\ni-j + y + 2x))^, g-\y) = y + {ln{j))^ . (5.14) 

Acting on solutions r, q or J, J of the Painleve IV equation with G raises u hj 2 
and keeps /i constant while g when applied on solutions expressed by j, J shifts both 
u and yU by 1 (see below) in such a way that acting twice with g agrees with the 
formula (I5.13p . Because of property of the g transformation to shift both parameters 
of the Painleve equation this transformation will be very useful in what follows in 
deriving solutions corresponding to new values of the parameters. 
Plugging r, q into expression (I5.10p leads to 

{lnq)^ = -j-j + j^/j, rq = ~p^ = -jj, (5.15) 
which after elimination of J yields a Riccati equation for j : 

= / + j(lng)x - . (5.16) 
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Rewriting equation (gZD as 2QP + + z/) = {2Pf + 2P(-Q - 2x) - 2P^ and 
recalling that Q = — (Ing)^. — 2x and 2QP + (/i + z/) = rq we find after comparing 
with relation fl5.16p that the solution to this relation is given by 

J = 2P = = ^ {y' -y. + 2xy - 2^) , (5.17) 

where use was made of the first equation in (14. 2p . Plugging the value of j from 
(I5.17P with two values of /i = ±|/i| into relations (I5.14p one obtains two set of trans- 
formations g± : z/^z/+l, /i^^(|/i|±l)^ and two set of inverse transformations 
gg^ : z/ ^ 1/ — 1, /i^ — >■ (|/i| ± 1)^ recovering expressions found in [8]. 

As observed in [1], if one sets Sn = (J) and i?„ = G" (j) then the above set 
of variables satisfy the equations of motion of the Toda chain: 

dx Sn = Rn+l ~ Rn^ Rn — Rn {Sn — Sn~l) ■ 

For quantities = J G" (J) dt and S^^lnr^ = J G" (J) dt one finds on basis of 
properties (15. 7p and (15. 8p of the G symmetry : 



G"(J) = Rn, (b.li 



from which the Hirota type of Toda chain equation: 



a^lnr. = ^^^±i^ (5.19) 



follows easily. The symmetry structure of g transformations is that of Volterra 
lattice [Ij. 

5.1 Darboux-Backlund transformations and "—2a; hierar- 
chy" 

Things simplify in the case of g = 0,r = r'^'^^ and thus with the initial AKNS Lax 
operator L^^^ = dx- A chain of Darboux-Backlund transformations define a string 
of new Lax operators via: 

= dx ^ L^'^ =dx + r^'^d-^^''^ 

with 



for k > 1. The eigenfunction r*^°) of has to satisfy 9j,^r*^°) = r*^°) = (9"r(°). 

The standard choice 

oo 
k=0 
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ensures that this is the case. To be able to assign a simple conformal weight to r 
one can break this sum and just set r^^\t) = Sn(t) due to the scaling law : 

5„(tA) = A"S„(t), tx = {\x,\%,...) 

which can be expressed as a differential condition (see ( 12.13^ ): 

d \ 

xd^ + 2t2-^ + . . . j Sn(t) = nSn(t), n = u-l. 

With this choice the Wronskian expression from [2] becomes : 

after using Sn-k = d^Sn- 
Similarly we find for g*^*^^ 

(jfc) _ 1 _ Wk-l[Sn-, Sn-1, ■ ■ ■, Sn^k+2] r ^ 

We see that the above Darboux-Backlund transformation changes the value of u by 
increasing the value of k. 

Note that for t2 = —1/4 and = = . . . = one gets 



00 



n.=0 

and so in this limit Sn{x) = Hn{x)/2"'n\. Thus r^^^ = Hn{x) and q^^^ = and 
1/ = n + 1 agrees with what we found above for solution with q = and u being a 
positive integer. 

Note that a simple Wronskian formula gives: 

j.{k)gik) _ ^k+l[Sn, Sn-1, ■ ■ Sn-k] Wk-l[Sn, Sn-1, ■ ■ •, Sn-k+2] 

W^[Sn, Sn-1, ■ ■ •! Sn-k+l] 
= —dl In Wk[Sn, Sn-1, ■ ■ ■, Sn-k+l] 

and so one finds 

p^"'^^ = 2d, In Wk[Sn, Sn-1, . . ., Sn-k+l] 

or in terms of the Hermite polynomials: 

p^'^'") = 2d, In Wk[Hn, Hn-l, . . ., Hn-k+l] , k>l, 

which satisfies p-equation fl2.16p with i' = n — 2-k + l and p^ = (ra + 1)^. That 
result is in agreement with 



(k,n) ^^ (k)\ n a l^ ^k+l[Hn, Hn-l, ■ ■ ■, Hn-kj , „ 

' > = (Inr'- ')^^-2x = d, \ In — ^ 1 - 2x 

y^ki-tin, J^n-1- ■ ■, J^n-k+l\ 
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satisfying (14. 4p with the same values of parameters. 

There exists another chain of successive DB transformations, which also builds 
the AKNS Lax pseudo-differential operator starting from the "pure" derivative 
by : 

with 

for k>l. The adjoint eigenfunction q'^'^^ of L*^°) has to satisfy dt^q^^^ = — (L^^^)^" g*^"-' 
(— l)"+^(9"g*^°\ The standard choice 



fc=0 



ensures that this is the case. To be able to assign a simple conformal weight for q it 
is convenient to break this sum by setting q^^\t) = Sn{t) due to the scaling law : 

With this choice the Wronskian expression for g'^'^^ becomes 

(fc) ^ Wk+i[q^^\d^q(^\...,d^,q(^'>] ^ W^fc+i 5^„_i, . . g„,fc] 
^ W^.[gW, 9.g(o), . . ., W,[S^, S^.„ . . ., 

after using Sn-k = d^Sn- Similarly we find for r^'''^ 

^(fc) ^ ^(fc) ^ 1 ^ iyfc-l[gn,'gn-l, ■ ■ ■,'?n-fc+2] ^ ^ 

Note that for t2 = —1/4 and = = . . . = one gets 



n=0 



and so in this limit Sn{x) = if„(a;)/2"n!, where Hn{x) is introduced in eq. (13. 4p 
Thus g^^-* = Hn{x) and r^^^ = and z/ = —n — 1 agrees with what we found above 
for solution with r = and z/ < and integer. Indeed 



p<'='") = 2d, In iy4^„, . . ., 



-fc+ij 



satisfies p-equation (I2.16P with z/ = — n + 2 ■ A; — 1 and /i^ = (ri + 1)^. 

To summarize we found the solutions to the p-equation (12.160 with u, /i param- 
eters given by : 

p^'^'") = 2d, In WklHr,., /f„_i, . . ., H^-k+i], u = n-2-k + l, fi' = {n + 1)' 

^^'^^ = 2d, In Wk[Hn, Hn-u . ■ H^-k+il u = -n + 2-k-l, fi^ = {n + 1)^ 

(5.20) 
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with n,k = 1,2,3, . . .. Corresponding solutions to the Painleve eqs. (14.31) and (14.41) . 
referred to as "— 2x-hierarchy" [5| [Bj are, respectively, 

^ik,n) ^ i^^(.-i) _ = 9. In ^,[H H „...,H ,^,] _ 2^ 

VVk-l[Hn, Hn-l, ■ ■ ., Hn-k+2\ (5.21) 

u = n-2-k + l, fi^ = {n + iy 

Wk-l[Hn,Hn-l, . . .,Hn-k+2] (5.22) 

V = -n + 2-k-l, ^^ = {n + Vf 
with n,k = 1, 2, 3, . . .. 

5.2 Darboux-Backlund transformations and "— 1/x-hierarchy" 

We will now examine an hierarchy of solutions which can be obtained by the 
Darboux-Backlund approach from remaining basic solutions listed in eqs. (13. 6 p 
and (13.51) . Namely, we consider as starting points q{x) = exp {—x'^) H^{x) with u 
being a positive integer and r{x) = e^^ H_^{x) for z/ < and integer. 
We begin with q^^^ = exp (— x^) Hm{x), r^^^ = and corresponding 



-d^ In (e ""^Hm^ - 2x = -d^ Inif^, 



which satisfies eq. (14.31) with v = m and /x^ = m?. 

The iteration procedure (15. 7p expressed here in terms of the AKNS variables : 



^« = G(g(°)) = -gW (lng(°))^^ + (g(°))V 



(0) 



(5.23) 



,(..1) = G(g(^)) = -g^'^) (lng(^))^^ + {q^^^Y ^ > ^ 

with g(°),r^°^ as defined above is solved by 

^(fc) ^ W^fc+ile""^ Hm-, e""" Hm+i-, • • Hm+k] 24) 

Using the recurrence relation (exp(— a;^)ifm)j; = — exp(— a;^)i7OT_|_i for Hermite's 
polynomials allows us to cast the above expression for q^'^^ into a simple form 

^ ^ ' Wk 

where 



Wk+i = Wk+i[f, fx, fxx,---Jx---x], f = e H, 



2 

m 
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Then applying the G transformation yields: 



(-1) 



k+1 



where use was repeatedly made of the identity WkWj! — (VT^)^ = Wk-iWk+i for 
an arbitrary argument /. Thus, G{q^''^) = 5^^+^) and (15.241) is established by an 
induction argument. 

Taking into account the identity 

iy,[e-^>i, . . ., e-^'Pfc] = W,[Pu . . ., Pk] e-^^' 

and setting m = n — fc + lwe obtain the following solution to the Painleve eq. (14. 3 p : 

v2 



Ak,n) ^ _Q Wk+l[e Hn+1, e ^ ifn, . . ., e ^ Hn-k+l] _ 2^ 

I? "2? T T 7" r _ T T -T' 2 7" 7" ^yi 2 t T" T 

(5.25) 



— 9 In ^'^1-'^"' -^"-1' ■ ■ ■' -^n-fc+i] 

with z/ = + + 1 and jj?' = {n — k + 1)^. Note that ?/(°'") reproduces solution 
—dx In ifji encountered above. 

We now turn our attention to r(x) = r'^'^\x) = Hm{x) and corresponding 

tf = c^a; In (^Hrri^ , which satisfies (14.41) with u = —m and /i^ = in agreement with 
eq. (13. 5p . Repeating steps applied above to find g'-'^-* one easily establishes that 



happens to be a solution to the iteration relations: 

r(-i) = G'-i(rW) = (inr^)^^ - (rW)^^^"), 
r(-^-^) = G-\r^-^y) = r(-^') (inr^''^))^^ + (r^-^))^ A: > . ^^'^^^ 

with q^^^ = 0. Plugging expression (I5.26P into expression for w from relation (14. ip 
and setting m = n — k + 1 one obtains solutions 

' ' = —Ox In '-^^ — = a^; In 



Wk+l[Hn+l, Hn, . . ., Hn-k+l] Wn^k+l[Hn, H^-l, . . ., Hk] 

(5.28) 
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to the Painleve eq. (14.41) with v = — (n + + 1), fi = {n — k + 1)^. In obtaining 
fl5.28p use was made of a Wronskian identity, 

Wk+l[Hn^l, Hn, . . ., ffn-fc+l] _ ^ ; -^n; ■ ■ -^fc+l] 

Wk[Hn,Hn-l,. . .,Hn-k+l] Wn-k+l[Hn, H^-l, . . Hk] 

where Cn,k are some combinatorial constants. 

Expressions fl5.25l) and (15.281) define "— l/x-hierarchy" [5]. 



5.3 Darboux-Backlund transformations and "—2a:/ 3- hierarchy" 
5.3.1 Solutions with fi"^ = (l/3)^ (2/3)2 

Consider now solution ( 13. 2p of the p-equation ( 12.16^ with /i^ = (2/3)^, = 0. Ac- 
cording to relation (14.111) the corresponding two solutions of equations (14.31) . (14. 4p 
are : 

y+ = (^x^ " -(Ing)x - 2x, y- = --^ Qa;^ ~ " (1^^)^ - 2x, 

with the following solutions of the string equations (I2.14p for u = 0: 

q = = --xe-2-'/3 = F^'^e-'^"/\ r = r^") = -xe^^"/'' = F^'^e^^"/' . (5.29) 
3 3 

In the above equation we employed notation involving polynomials : 

defined for = 0^ 1, 2, n = 0, 1, 2, . . .. By simple rescaling of their arguments poly- 
nomials F^\Fn become proportional to Hermite polynomials Hm-,F[m for certain 
values of m. 

Acting on the initial configuration (I5.29p successively with the DB transforma- 
tions as in (I5.23P one arrives at the following Wronskian representation in terms of 
the ratios of Wronskians : 

2 

/i = ±-, u = 2n, n>0. 
3 

Acting with negative powers of G as in (I5.27P yields : 

W 9(1) 



wM'\F['\...^Fi 



^^exp(|!], g(-) = -l/r(-'^+^) (5.31b) 



2 

ji = ±-, v = —2n, n >0 . 
3 
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The first few solutions are 

^^ 2(-12x2 + 4x^-9) _2x! 8x(504x^ - 192x^ + IGx^ - 2835) _2xi 

' = — 3 q^^) = ^ ^ ie 3 

^ 27x 243(-12a;2 + Ax'^ - 9) 

t ^^ 2(12x2- 9 + 4a;^) 2^ 8x (-2835 + 504x^ + 192 + 16 a;^) 

r*-^ = — ^ -e 3 ^i-^J = 1 ig 3 

27x ' 243(12x2 - 9 + 4x4) 

Applying recursively (15. 3p one gets closed Wronskian expressions for the p- 
f unction: 

4x 

y 

4x 



pin) ^ Qu^p^ =p-2n^ + 2 (in W^[Fi,'\Fi'\ . . F^i^ 
= ^-"(p) = p + 2n^ + 2 (in . . ., 



which are solutions of the p-equation (12.161) with u = ±2n, n = 1,2,3,..., respec- 
tively. In particular 

p(±i) = J_(4a;4_^27±36x2) 
27x 

satisfies the p-equation (I2.16P with p"^ = |, = ±2. 

Consider now the basic solutions (13.11) of p-equation (12.161) . Plugging (13. ip (with 
a plus sign) into relation (14. lip one obtains the following two solutions of equation 
(SD: 

2 1 , , 

y = y+ = -^x, ^ " ^3' ^ = l^^-^^) 

W = y- = -7T^ 16 o , 8 ' /^=±7T, v = ~l (5.33) 

O 3 



Eq. (I5.33P corresponds to the following AKNS variables: 

3 V3 ; ° (5.34) 

q = = _e-2-'V3 ^ _i.(o)e-2-'V3^ 

Applying the DB transformations G, generalizes the above basic solutions to 



g^"^ = - ^"^^^!^^, '2 '""m ^ ^^P(-^K r("-^^) = l/g(") (5.35a) 
W4Fl,'\Fi'\...,Fi%] 



1 

ji = ±-, = — 1 + 2r;,, > 
3 



and 

ri-n) 



W4Fi'\Fi'\...,F^'^ 



'] /2t2\ 

-i exp ^ , g(-") = -l/r(^"+^) (5.35b) 
1] V 3 / 



1 

/i = ±-, v = —1 — 2n, n > . 
3 
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Moreover one obtains the following Wronskian expressions for the chain of associated 
solutions to the p-equation: 

pin) = G'"(p) =p-2n^ + 2 (in W^[Fi'\ F^, . . ., Ff j)^ , u=l + 2n, 
pi-n) ^ ^-n(^) = p + 2n^ + 2 (in W^F^'Kf^, . . ., F^Jj)^ , u = 1 - 2n 
with n = 1, 2, 3, . . .. 

Plugging (13. ip (with a minus sign) into relation (14. lip one obtains the following 
two solutions of equation (I4.3P : 

2 3 + 2a;2 1 ^ ^ 

y = y^ = -r^3T^^^ /^ = ±3, -=+1 (5.36) 

2 1 , , 

w = y- = --x, '""^3' ^^=+1 (^•2'^) 

which corresponds to the following AKNS variables: 

r = r(°) = F(i°)e2^-'/3, q = q(^^ = -F^'^ e-^^"/\ /i = ±-,j/=l. 

3 

Applying the DB transformations G, G^^ generalizes the above basic solutions to 

„W _ Wn+l[Fl^'^\Fi^\...,Fi^^] f_'^\ (n) _ . I (n-l) /p- oo„X 

/i = ±-, z/ = 1 + 2n, n > 
3 



and 

^(0) f(0) p(0)i /o^2 

-n) _ 

1 

/i = dz-, z/ = 1 — 2n, n > , 
3 

The Wronskian identities : 

W^,,[f!^\ Fr\ . . ., = -l^„[i^i^ F?\ ■ ■ i^S], 

w^,AFis'\ f['\ . . ., m = wM'\ Fr\ . . 

ensure that the solutions w of the Painleve IV equations generated by expressions 
in eqs. (I5.38P coincide with those solutions w which originate from equations 
(I5.35P for equal parameter v (accomplished by shifting n to n ± 1 when going from 
(I535|) to dSSH])). 

Moreover one obtains the following Wronskian expressions for the p-function: 

p'-' = G-{p) = p + 2nf + 2 (in W^[Fi'\ , • • ^fj 

Ax 



, n = 1,2,3,... 



p(«) = G"(p) = p - 2n^ + 2 (in W^F^'K F^ , . . ., Ff_V 



n= 1,2,3, 



which satisfies the p-equation (12.160 with p^ = |, = 1 ± 2n. 

Wronskians Wn[Fi'^\Fi^\ . . ., F^^^] and Wn[F^°\ F^^^K . . ., F^^^] are proportional 
to the Okamoto polynomials Qn as defined in |5j. 
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5.3.2 Solutions with fi = ±4/3, ±5/3, . . . 

To reach expressions for solutions with n = ±4/3, ±5/3, ... we have employed trans- 
formations g, with their properties of raising and lowering fi and u by one when 
applied to j,j configuration. This property of the g,g~^ transformations ensures 
that we reach all allowed values of the Painleve IV parameters following the zig-zag 
DB orbits through the (/x, z/)-plane. 

We now present results obtained by applying g, g"^ transformations to the three 
basic cases with the values /i^ = (2/3)^, (1/3)^ presented in equations (15.311) . fl5.35p 
and (15.381) in the previous subsection. 

We start with a chain of solutions (I5.3ip to the Painleve IV equations with 
/i^ = (2/3)^. Through the successive actions of g, g~^ transformations these solutions 
generalize to : 



fi^ = [ — |-?n±l) , u = 2k — m — 1 





(5.39a) 



2 



f/=(- + m + l\ , u = l + m-2k 




and solve the Painleve IV equations (14. 3 p and (14.40 . respectively, for the positive 
integers k,m > 0. 

Generalizing through the DB approach the system of solutions given in eqs. 
(I5.35P leads to : 

u/ rp(i) p(o) p(o) p(o)i 

_ I In I '^fe+"^+2[-^o ^-^1 ^•••^-^"^ ^-^0 ^-^1 ^■■■^-^k \ 1 1 _ n (K AC).] 
' ' M/ IpW p(i) p(o) p(o) p(o) ill o ' lo-^uaj 




2,1 \^ /2 , 

/i = (— -±m±lj = l- ±m) , u = 2k — m — 2 



/i^ = l^-i ±m± 1^ =Q + m^ , z/ = -2A;±m±2 

written in terms of two positive integers k and m. 

By applying the same procedure to solutions given in equations (15.381) one obtains 
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a new class of solutions 



= [ — \- m ] , u = 2k — m — 1 



2 




.2 ' 1 



/i = + mj , u = m — 2k + 1 

written in terms of two positive integers k,m>0. 

It is important to point out that the first two of the above three classes of 
solutions overlap in their values of the fi and u parameters for /x^ > (2/3)^ and 
in such cases the corresponding solutions coincide. For instance, solutions fl5.39ap 
with /i^ = (2/3 + m + 1)^, u = 2k — m — 1 are equal to solutions fl5.40ap with 
/i^ = (—1/3 + m' + 1)^, u = 2k' — m' — 2 when m' = m + 1 and k' = k + 1. 
Thus, for any solution ym,k or Wm,k from f l5.39p there exists an identical solution 
ym'=m+i,k'=k+i OY Wm'=m+i,k'=k+i from fl5.40p . As cucountered before in the text, an 
equality between two classes of Painleve IV solutions derived by action of the DB 
transformations stems from existence of a special Wronskian identity. In this case 
the relevant identity is given by : 

i^.+..2 F«, . . ., , f( vr , • • i^r] = 

— Kk^myVk+m+AyrQ , . . ., i^o i^i , • • -f^fc+ij , 

with certain combinatorial constants Kkm- 



6 Conclusions and Outlook 

We presented here a systematic and self-contained derivation of rational solutions 
to the Painleve IV equation using the method of the Darboux-Backlund transfor- 
mations of the particular reduction of the AKNS pseudo-differential Lax hierarchy. 
By studying the orbits of the Darboux-Backlund transformations originating from 
few seeds solutions we were able to find closed expressions for solutions associated 
to all allowed values of the Painleve IV parameters. The explicit expressions for 
all Wronskian representations derived here seem to be in general agreement with 
determinant formulas obtained earlier in [TOl [121 E] by alternate methods. 

In a separate publication we will discuss the Darboux-Backlund transformations 
versus the (affine) Weyl group symmetry of the Painleve IV equation [13j . 
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